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Abstract 

We present the results of studying the influence of different environmental 
states on the coherence of quantum processes. We choose to discuss a simple 
model which describe two electronic reservoirs connected through tunneling 
via a resonant state. The model could, e.g., serve as an idealization of in- 
elastic resonant tunneling through a double barrier structure. We develop 
Schwinger's closed time path formulation of non-equilibrium quantum statis- 
tical mechanics, and show that the influence of the environment on a coherent 
quantum process can be described by the value of a generating functional at 
a specific force value, thereby allowing for a unified discussion of destruction 
of phase coherence by various environmental states: thermal state, classical 
noise, time dependent classical field, and a coherent state. The model allows 
an extensive discussion of the influence of dissipation on the coherent quan- 
tum process, and expressions for the transmission coefficient are obtained in 
the possible limits. 
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I. INTRODUCTION 



The progress in fabrication of submicron structures has lead to a wealth of new structures 
whose transport properties are dominated by the feature that electronic transport through 
the structure takes place coherently. However, additional degrees of freedom in addition to 
those of single electrons are present and it is of importance to account for their influence on 
coherent processes. In the following we shall investigate how different environmental states 
influence the quantum process of main interest. Although we shall choose a model which is 
directly relevant to transport properties of a double barrier structure or transport through 
a quantum dot, the technique we develop is of general interest. 

The paper is organized as follows: In section || we model the system whose transport 



properties will be studied, and in section 111 we show how to describe the transmission 



probability of a structure in terms of Green's functions. In section |V| we discuss the main 
approximation, which will allow an analytical description of the influence of dissipation on 
the transmission properties. In section we introduce the closed time path formulation 
and develop the description of the destruction of phase coherence in terms of a generating 



functional technique, and in section |V| we discuss the effects of different environmental 
states: thermal and coherent states, the external field case, and the fluctuating level model. 
We shall for all these cases derive closed expressions for the transmission probabilities, and 
compare the results in order to notice similarities and differences in the influence of different 
environments. In section |VII| we discuss destruction of phase coherence using the generic 



Aharonov-Bohm situation. Finally in section [VIII| we summarize and conclude. 



II. THE MODEL FOR STUDYING DESTRUCTION OF PHASE COHERENCE 

In this section we shall set up a minimal model describing the dissipative feature of 
the general quantum transport problem. The simplifications are introduced in order to be 
able to treat analytically the influence of environments on the transport process. We shall 
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study the general problem of quantum transport where electronic current reservoirs provide 
electrons to an active region where interaction with an environment can take place. The 
reservoirs corresponds in reality to large electrodes, and we can describe the Hamiltonians 
for the left and right electrodes in terms of their electron energy levels 

= 5Iep,aaJ,_^ap,^. (1) 
p 

The quantum number p labels the momentum of the electron states, and a = l,r refers 
to the left and right electrode, respectively, and aj^^^ creates particles in states with these 
quantum numbers with the corresponding energy ep^o-. 

The central active (sample) region can in the absence of coupling to the reservoirs be 
described by a Hamiltonian or equivalently its energy levels. The single particle energy at 
the central site, labeled c, is ec corresponding to the term in the Hamiltonian 

Hs = ecalac- (2) 

In the following we shall consider a model which allows extensive analytical calculations, 
and therefore restrict the number of levels relevant in the central region to one. Eventually 
we shall consider also the case of two levels. In the event that the reservoirs are connected 
through the central site, electrons are transmitted between them. Such a situation can 
quite generally be modeled by transfer matrix elements, V^^o-, between the reservoirs and the 
central region. The coupling of the electrodes to the central site is therefore described by 
the tunneling HamiltonianS 

= E{'^P'x4«P<x + h-c.}. (3) 

PCT 

The transfer matrix elements are here considered to be phenomenological parameters, but 
can of course for any chosen microscopic model of say a double barrier, be expressed in 
terms of the potential profile and the carrier mass. The Hamiltonian, He, for the electronic 
system of interest is therefore given by 

H, = Hi + Hr + Hs + Ht. (4) 



Within the sample region we allow for interaction with an environment. For our purpose 
we can quite generally assume a bosonic environment with a corresponding Hamiltonian Hj, 
which has the standard normal mode form 

H, = Y,huj^i^bib^ + ^'j (5) 

in terms of the bosonic creation and annihilation operators, b^^ and ba- 

For the coupling of an electron at the central site we take a linear coupling to the normal 
modes 

Hi = alttcX, (6) 
where X is the collective environment displacement operator 

X = J2X.{bi + b^}, (7) 

a 

and Xa the coupling constant to mode a. 

The resulting Hamiltonian for system and environment 

H = H, + H, + Hi (8) 

has been discussed in a variety of contexts!, most recently in the context of inelastic reso- 
nant tunnelingH, where the model has served as a simplified description of the influence of 
interaction with phonons on the transport properties of a double barrier structure. In the 
present account we shall not only discuss a thermal environment, but a variety of environ- 
mental states and their influence. A purpose of the paper is to present a calculational scheme 
that allow a unified description of arbitrary environments, and in a direct way exhibits the 
physical content of a model so that calculations within more realistic models can be made 
tractable. We now proceed to describe in detail the transmission properties of the system 
under consideration. 
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III. TRANSMISSION PROPERTIES 



In mesoscopic physics where the main feature of electronic transport is its coherence, the 
transport description can be expressed in terms of the scattering properties of the mesoscopic 
structure. The choice of model has been dictated by this feature, and in the following 
we shall study the quantum mechanical problem of transmission of an electron through 
a region where it can interact with additional degrees of freedom. We wish, therefore, 
to calculate the transmission coefficient for a particle emitted, say, from the left reservoir 
under the assumption that it propagates via the central site, where it is allowed to interact 
with an environment, to the right reservoir. In accordance with the initial condition of an 
electron impinging from the left reservoir, we can assume that at some initial time ti the two 
subsystems, particle and environment, are decoupled, so that the initial state is described 
by a separable statistical operator 

Pi = Pp',i * Pb, (9) 

where Pp'^i is the projection operator describing an electron in the left electrode in momen- 
tum state p' 

^'pv = Ip',0(p',^I- (10) 

The state of the environment is presently arbitrary and described by the statistical operator 
Pb. As regards the transmission problem, given the above initial state, all information can 
be extracted from the probability, Pp' ^i^p^r{t) , to find the electron in the right electrode at 
time t with momentum p. This conditional probability is given by the expression 

Pp',l^pAt) = MP^UHt, ti)Pp,rU{t, U)), (11) 

where Pp^r is the projection operator 

Pp,,= |p,r)(p,r| (12) 

corresponding to the assumed final outgoing particle state with momentum p in the right 
electrode, and 
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C/(t,t.) =e-t^(*-**) (13) 

is the evolution operator corresponding to the total Hamiltonian. The trace with respect to 
all the degrees of freedom is denoted by Tr. The absence of any environment operator dis- 
criminating the final states of the environment is in accordance with the typical experimental 
condition pertaining to the electronic conduction process, namely that the environmental 
degrees of freedom are not observed. 

For an electron to propagate between the reservoirs it first has to enter the sample region 
and lastly to exit it. In the event we only explicitly consider interaction with the environment 
in the sample region, we can exploit this feature and introduce a discussion in terms of the 
effective coupling between reservoirs and sample region. The chosen model has this feature 
to the extreme. In the assumed model the electron only couples to the environment at the 
central site. We can therefore express the transition probability in a form that explicitly 
only involves the dynamics of the electron at the central site and the environment by simply 
noting, that in order to calculate the amplitude for a transition from the left to the right 
reservoir, the first propagation has to be from the left reservoir to the central site, and the 
last propagation from the central site to the right reservoir. The transition probability can 
therefore be rewritten as 

Pp',i^p,r{t) = Al^riVd" I'dti fdt2 fdh fdU 
h Jo Jo Jo Jo 

exp |^ep,r(t2 - h) + ^ep/,/(t4 - 

{Gf{U,ts)G^{h,t,)). (14) 

where (...) = tr{pb...) is shorthand for the trace with respect to the environmental degrees 
of freedom, weighted with respect to the initial environment state, and 

G^{t,t') = -i9{t-f){0\[a,it),alit'm (15) 

is the retarded Green's operator (operator with respect to the environmental degrees of 
freedom) for the central site dynamics as 
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a,(t) = e^^*a,e-^^*, (16) 

and |0) denotes the particle vacuum state. The advanced Green's operator is given by 
hermitian conjugation 

Gf(t,t') = [Gf(t^t)]^ (17) 



where f denotes hermitian conjugation. In eq.(H), the arbitrary initial time has been chosen 
at time zero, ti = 0. 

In many cases of interest we do not need the full information on the transition probability 
as function of time. For instance, if we are only interested in the average particle flow in a 
steady state we are only interested in the transition probability per unit time 

Wp,,^p, = hm ^El±:^^ (18) 

or the transmission coefficient for making a transition between a state of energy e' in the 
left reservoir and a state of energy e in the right reservoir 

T(e, e') = h^Wpi^i-,p^r5{e- ep^r)Kt - tpn). (19) 
pp' 

Invoking the scattering approach for the description of transport properties for coherent 
quantum processes! we have at zero temperature the contribution to the conductance at 
energy e' 

G(e') = ^T(e'), (20) 

where 

T{e')= / d6T(e,e'), (21) 

<^ 

is the total transmission coefficient, the probability to reach the right reservoir for a state of 
energy e' in the left electrode. The factor of two is the spin degeneracy factor of the electron. 
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IV. THE WIDE BAND APPROXIMATION 



We have in the previous section reduced the expression for the transition probabihty to 
an expression which only involves the dynamics at the central site. However, since both 
tunneling and interaction with the environment is present the dynamics is complicated 
and no closed expression for the transmission probability can be found, that is, without 
any further assumption no simplification is possible as the phonon average cannot be done 
explicitly. If we, however, assume that the widths, W, of the electronic energy bands of the 
electrodes are the largest energy in the problem, W > TitOa, it is possible to obtain closed 
expressions for the transition probabihty for various cases of environmental states. In this 
wide band limit an electron at the central site decays into a continuum of states in the 
electrodes, and we expect to have a quasi-stationary state which decays exponentially in 
time. To utilize this property of the electrodes we first discuss the lifetime for occupation of 
the central site in the absence of an environment. 

In the absence of coupling to the environment the retarded Green's operator for the 
central site, G^, reduces to the c-number Green's function, equaling the amplitude for 
the particle to remain at the central site. The Fourier transformed amplitude satisfies the 
Dyson equation 

G^{e) = g^{e) + (6)Ef (e)^^ (e), (22) 
where the central site self-energy is given by 

^^{^) = Y.KA'9U'l (23) 

and describes the coupling of the central site to the reservoirs, as indicated by the appearance 
of the electrode propagators. A lower capital is used to designate a subsystem Green's 
function, implying that the electron is not allowed to propagate in or out of the electrodes, 
i.e., its dynamics is determined by the Hamiltonian with all connecting elements, V^,(7, set 
equal to zero. 
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In the absence of tunneling between the electrodes and the central site we have for the 
isolated electrode Green's functions 

and for the central site Green's function 

In the presence of coupling between electrodes and the central site the Dyson equation yields 
the central site Green's function in terms of the self-energy 

The escape rate from the central site is given by the imaginary part of the self energy 

r(e) = -$>mEf (e) 

-7r^|l^p<,|25(e-ep,). (27) 

If the hopping matrix elements vary slowly with energy in the resonance region, e ~ Sc, 
and the energy bands in the reservoirs are wide so that also the electrode density of states, 
N^{e) — J2p ^i^~^pa), has a weak energy dependence, we can neglect the energy dependence 
of the escape rate, and by analyticity we have 

- Loo V e-e' + iO' ^^^^ 
which implies that the real part of the self energy vanishes in the wide band limit. 

The wide band approximation simplifies the problem considerably as the time dependent 
Green's function becomes a decaying exponential 

Gf{t, t') = -ie{t - i')e-t(^c-ir)(t-t')^ (29) 
and the Green's function satisfies, for times t > t" > t' , the relation 

Gf{t,t")G^{t",t')^G^{t,t'). (30) 



This group property leads to a tremendous simplification of the interacting problem as the 
tunneling dynamics at the central site and the environment dynamics decouple. If the group 
property, eq.(pOD, is valid we have for the central site amplitude correlations 

(G^(t4,t3)Gf(t2,tl)) 

= G^{t4,t3)Gj^{t2,ti)Z{t4,ts,t2,ti), (31) 

where we have introduced the influence function describing the effect of the environment 

ZiU.t.MM) = {{fe^^n"' ^W)(re-^^''* ''^'^)) (32) 

with the environment variable in the interaction picture 

X{t) = e^^^*X e-^^"*. (33) 

and T and T denotes the time and anti-time ordering operators, respectively. The decoupling 
of the particle and environment degrees of freedom is visualized using Feynman diagrams in 
figure 0. 

Because of the group property, eq.(pOD, the central site dynamics has a special behavior: 
the lifetime of the central site state is independent of the coupling to the environment, as 
the probability, Pc{t), for the particle to remain at the central site after a time span t, if 
initially at the central site, is 

Pe(t) = (Gf(0,t)Gf(t,0)), (34) 
and noting that Te~ft -^o '^^^(^ jg the interaction picture time evolution operator for the bath. 



and consequently unitary, we have Z(0, t, t, 0) = 1. Therefore by eq. (|3lD we have for the 
staying probability 

P,(t) = (Gf(0,t)Gf(t,0)) = e-f*, (35) 

which is independent of the coupling to the environment .@ 

Before we go on to calculate the influence function Z, which contains all information 
about the influence of the environment for various environmental states, we briefly discuss 
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the transmission problem in the absence of couphng to the environment. We shall explicitly 
assume the exponential form eq.(p9D for the decay amplitude. 

If we for the moment neglect the coupling to the environment we have for the transition 
probability 



i|Vp,rr|Vp',,|' f'dt, fdt2 fdh f'dU 
n Jo Jti Jo Jo 



g|ep,r(i2-t3) + ^ep,i(*4-tl) 



e 



-U^c-ir){t2-ti)+j-{e,+ir){t3-u)_ ^3g-j 



The integrations are readily done and we obtain for the transition probability per unit time 



w^p',/^p,r = lim 

IT/ |2|T/,,|2 

S{^p,r - ^p',l), (37) 



t^oo f 

27T |Vp,,|2|Vp,z|2 



or equivalently for the transmission coefficient 



9F r 

^'''^A,(e')5(e-e'), (38) 



r 

where we have introduced the central site spectral weight function 



2V 



and the left and right escape rates 



(39) 



= vr^|Vp,p5(e-ep.). (40) 



p 



We observe that the transmission coefficient has the expected resonant character of the 
Breit-Wigner formula^. 

In the absence of coupling to the environment we therefore have the contribution to the 
conductance at energy e' 
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G(e') 



4e^ r,r, 

T:h (e'-ec)2 + r2' 



(41) 



Having discussed briefly the uncoupled case, we now turn to calculate the influence 
function Z for various environmental states. We shall invoke the assumption of a wide band 
width allowing for the decoupling expressed by eq.f^l]), and obtain in this approximation 
the following expression for the transmission coefficient 



In the wide band limit we notice the integral of the transmission coefficient is related to the 
staying probaility at the central site 



and we can therefore conclude that the integral of T(e') is unaffected by the coupling to the 
environment since the staying probability, Pc{t), is unaffected. 

V. CLOSED TIME PATH FORMULATION FOR THE INFLUENCE FUNCTION 

In order to calculate the influence of the environment on the otherwise coherent quantum 
process, the task is to calculate the influence function given in eq. (p2D . The flrst step in 
achieving such a goal is to transform the operator expression into an analytical one. A 
general and convenient method for performing this transformation has been devised by 
Schwingeiil. The key point to note for employing this method for the present purpose is to 
note that the influence function Z quite generally can be expressed in terms of a closed time 
path generating functional provided the "force" is chosen properly as 






(43) 



(44) 
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where C is the closed time path extending from — oo to +00 along Ci and back again along 
C2, as illustrated in figure ^ and Tq is the contour ordering symbol, ordering the environment 
operators X{t) according to their contour label position on the contour C (earliest to the 
right) and /°(t) is the function on the upper and lower branches of the contour specified by 



is therefore the quantity of interest. 

The infiuence of the environment only appears through the state of the environment in 
the generating functional, and different environments can now be handled on an equal footing 
corresponding to just substituting the proper environmental state, i.e., the corresponding 
statistical operator. 

In the present formulation we note that corresponding to the infiuence of an environment 
on a given physical quantity there corresponds a function, here denoted the proper "force", 
which when inserted into the generating functional completely specifies the infiuence of the 
environment. 

In the following section we calculate the generating functional and the infiuence function 
for various environmental states. 

VI. THE INFLUENCE FUNCTION FOR TYPICAL ENVIRONMENTAL STATES 

In the preceding section we showed that in the presented model the infiuence functional 
refers only explicitly to the environmental degrees of freedom and is characterized by a single 
"force". This simplifying feature will allow us to obtain closed expressions for the infiuence 
function for typical environmental states. 




(45) 



The closed time path generating functional 



z[f] = {Tce-Uj-n-)nr)^ 



(46) 
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A. Thermal environment 



We first discuss the quite common physical situation where the environment is acting 
as a heat bath. This could, for example, be the case in question where the Hamiltonian is 
thought to represent a resonant tunneling structure, the lattice degrees of freedom of the 
crystal acting as the heat bath. 

A heat bath is characterized by a single macroscopic parameter, its temperature T, and 
the environmental state is in this case described by the equilibrium statistical operator 

P^ = PT = tre-H,/k,T - (47) 
The average of the generating functional is then Gaussian, yielding the quadratic form 

Z[f] = exp[-^ Jdr Jdr' f{r)D{r,r')f{r')Y (48) 

where 

D{T,r') = -z{Tc{X{r)X{r'))) (49) 

is the contour ordered bath Green's function. 

It is convenient for the physical interpretation to split the exponent appearing in the 
generating functional into real and imaginary partJli 

Z[f] = exp - -2 / dt dt' [2f^{t)D%t, t') Uit') 

L ft, J— CO J —oo 

+/_(t)D^(t,0/-(i')]}, (50) 



where 



f±{t) = \{h{t)±mi (51) 



and 



D\t,t') = -ie{t -t'){[X{t),X{t')]) (52) 
is the retarded bath propagator, and 
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D^{t,t') = -z{{X{t),X{t')}) (53) 

is the correlation or Keldysh bath propagatoi0. 

For the present thermal case we have for the retarded Green's function 



h roc 

D^{t,t') = 0{t-t') du;J{u;)smuj{t-t'), (54) 

2 Jo 



where 



Ji^) = ^ E >^'Ji^ - = -^'^mD'^iuj) (55) 

is the spectral function completely characterizing the influence of the microscopic degrees of 
freedom of the bath. For the correlation function we similarly have 

Ti^ Tloj 
D^{t,t') = -i— duj J (lu) coth—— cos uj{t-t') . (56) 
2 Jo zkgT 

The Fourier transforms of the retarded and the correlation propagator is connected according 
to the fluctuation-dissipation theorem 

D^(cu) = 2z53mD^(cu) coth (57) 

The calculation of the influence of the thermal reservoir on the transition probability 
has now been reduced to the performance of simple integrals which are readily done and we 
obtain for the influence function 

Z (t4, ^3, ^2, ^l) = Z^(t4, ^3, t2, ti) Z^(t4, ^3, ^2, ^l), (58) 

where the contribution from the retarded bath propagator is given by 

Z«(t4, ts, t2, ti) = e3^(*=^-*^+*^-*3) Z^(t4, h, t2, tl), (59) 

with 

Z'^iU, t2, h) = exp <^ - dio ^S^iU, h, t2, ti) \ , (60) 
specified by the function 
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Sui{t4:t3,t2,ti) = —smu;{t2 — h) — smu;{t3 — ^2) 
+ sma;(t4 — ^2) + sinci;(t3 — ti) 

— sina;(t4 — ii) + sina;(i3 — ^4), (61) 
and the effective coupling constant given by 

A=rd.M. (62) 

Jo LO 

The contribution from the correlation bath propagator is similarly given by 

rvK , \ f 1 /■°°, , 

Z {t4,t3,t2,ti) = exp\ - - duj — 5-coth— — 

I 4 JO oj"^ zkb-i 

Co,ih,h,t2,h)y (63) 

where 

Cu,{t4,t3, t2,ti) = 2 - cosa;(ii - ^2) - cosa;(t2 - ^3) 
+ cosa;(ti — ts) + cosa;(t2 — ^4) 

— cosa;(ii — ^4) — cosa;(i4 — ^3). (64) 
The two real functions and S^^ are the real and imaginary parts of the complex function 

UU, h, t2, ti)^2- e*-(*2-ti) _ ^Mt3-t,) 

_|_gia;(t3-ti) _|_ gia;(t4-t2) 

Since the environment is in thermal equilibrium the influence function ^^(^4, ^3, ^2, ^i) in 
fact depends only on three independent variables 

r ^t2-ti 
t' = h- U 

r"^h-t2, (66) 

the above choice being determined by the original time labeling of the Green's functions. 
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In terms of the independent variables we then obtain for the transition probabihty per 
unit time 

y Pit) 
\*'p,r\ rp',Z| / 7 _// / 7 / / 



dr" / dr' / dTZ{r,r',r") 

-OO Jo Jo 

exp {-^(w - - ^(v,; - - ^r)'^' 

+^(epV-e; + ir)T|, (67) 

where the influence function in terms of the three independent variables again is split into 
the two distinct parts 

Z{r, r', t") = Z^(r, r', r") Z^{t, r', r"), (68) 



where 



with 



Z«(r,r',r") =exp|^_^°^ciu;^5.(r,r',r")}, (69) 

Suj{t, t\ t") = — sina;T + sina;T' 

— sinter" + sina;(T" — r') 

+ sina;(r" + r) - sina;(r" + r - r'), (70) 



and the correlation part 



Z^(r,r',r") =exp|-i Ho. ^C.(r, r', r") 



where 



coth^l, (71) 
2kBT J ' ^ ^ 



C^{t, t', t") = 2 — cosa;T — cosa;T' 

— cosa;T" + cosa;(T" — r') 

+ cosa;(T + t") - cos uj{t" + r-r'). (72) 
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The central site energy is shifted downwards according to 

ic = ec-^ (73) 

similar to the negative polaronic energy shift. 

For the transmission coefficient in the case of a thermal environment we then get 

9 /"OO fOO POO _ 

T(e,e') = ^r^r, / dr" / dr' / dTZ{T,T\T") 
irh J-oo Jo Jo 



exp 



+ -(e'-e~ +zr)r|. (74) 

This expression has been studied perturbatively in the case of the Einstein mo dei, and also 
in terms of elastic and inelastic channels^. In the following we investigate the total trans- 
mission coefficient in the thermal case and obtain explicit expressions for various limiting 
situations. 



1. The total transmission coefficient 

If we are not interested in the energetics of the arriving particles in the right electrode, 
but only in the number of arriving particles, only the total transmission coefficient is relevant. 
Such a situation arises for instance in the case where we can neglect any effect of the Pauli 
principle in the right electrode, corresponding to the situation of a highly biased left electrode 
in which case the left-going current is zero. In this situation the current of arriving particles 



IS 



roo 

I = e de'T{e')hie'). (75) 
Jo 

In such a case we only need the expression for the total transmission coefficient which. 



according to eq. (|7^ , pID , is given by 

4 



T(e') = T^r^r, / dr' drZ{T,r') 
n Jo Jo 

exp |-^(e' - e, - zF)/ + ^(e' - + zr)r| 



(76) 



where 



Z(r,r') = ^^(r,r')^^(r,r'), 



(77) 



is specified by 



Z«(r,r')=^^(r,r',0) 



exp 



A Jo 



(78) 



and tiie correlation part 



Z^(r,r') = Z^(r,r',0) 



~, J(a;) , Tiuj 
\ cotn — — — 

1 — cosa;(r — t'))|. 



(79) 



Since the environment is in the thermal equilibrium state, the influence function for the 
total transmission coefficient only depends on one time variable, and the expression for the 
total transmission coefficient can be reduced to a single integral by introducing the mean 
and relative time variables 



r + r' 



(80) 



For the integration region we observe 



POO POO roo roo 

Jo ^^io ^^' = loo^'4 



(81) 



and performing the integration over the mean time we finally obtain 



Z^{t) exp 1 ^(e' - Qt \ + cc. 



(82) 



where 
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Z^{t) = exp <^ - 



4 Jo 



~, J{u) . 
duj — — sm ut 



and 



Z'^W=exp|-l/ 



1 /•°°, Jiuo) ^ huo . ,1 
duj — — cotn — — — (1 — cos cut) > . 



(83) 



(84) 



4 7o 2kBT 
In the case where the oscillators all have the same frequency, the Einstein model, the 
spectral function takes the form 



4A^ 

n 



(85) 



relevant, e.g., to optical phonons. In this case, we deduce from eqs. ( ^^ - |8^) the total trans- 
mission coefficient 



T(e') 



2TlTr 



exp 



dt exp 



An 



An 



hujQ 



l+n(cuo))e-*"»* 



1 + 2n{uJo)) 



exp 



'2^ 

flUJn 



n{ujQ)e 



iujQt 



exp |-^(r - i{e - ec))t| + c.c. 



where n(uj) is the Bose function 



n(uj) 



1 



(86) 



(87) 



Expanding the exponential functions containing e*"^"* and e"*"^"*, and performing the inte- 
gration over t we obtain 

Pn{T) 



(e' - e, - nhuj^Y + 
where the temperature dependence is specified by 



(88) 



Pn{T) = exp 



E 

ni,n2 = 



V 



:0 ni\n2\ 



[1 + 2n(^o)] 



huJn 



[1 + n{uoo)\ 



n2 



n{uJo) 



(89) 
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the prime restricting the summation to the terms for which rii — n2 = n. Transmission 
can take place with absorption or emission of oscillator quanta giving rise to additional 
resonance peaks, the function Pn{T) determining the relative weight of the peaks. The 
energy dependence of the transmission coefficient is illustrated in figure 0. In particular, in 
the low temperature limit, ksT < huQ, where the Bose function vanishes, n{uj) 0, we can 
ignore all terms containing powers of n{uj) in eq.(|89|), and obtain at zero temperature 



exp 



\hujoJ f n! \nujoJ 



P„(0) = <( I ^""^"^ J " ^""^"^ , (90) 

n < 

reflecting the possibility of an electron off resonance to tunnel from the left to the right 
reservoir via the central site by spontaneously emitting a phonon, and the impossibility 
of gaining energy from a zero temperature environment. We note that -P„(0) is a Poisson 
distribution characterized by its mean value 



(n) = f;nP„(0)=fA) . (91) 



n=0 



2. Approximations conserving the integrated transmission probability 

To illustrate the general features of the systematic and fluctuation influences of the 
environment we choose the spectral function in the further calculations to have the form 

JM = ,.(^)'"exp(-5^), (92) 

We note, that rj is a dimensionless constant describing the coupling strength between the 
central level and the environment, and Uc is the cut-off frequency for the oscillators. 

There are several limits in which simple expressions for the transmission coefficient may 
be worked out. These limits are the broad resonance limit, F > huc, the strong coupling 
limit, r] > 1, the high temperature limit, ksT > (1 + l/ri)huJc, and the weak coupling and 
low temperature limit, i] < 1, ksT < — . 
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a. Broad resonance limit. In the broad resonance limit, F > TiuOc, a short time expansion 
of the influence function is sufficient. The correlation part of the influence function is 
therefore given by 



(t) = exp 



which controls the other part to become 



(93) 



A 



Z«(t) =exp -^-t . 



(94) 



Here k is defined as 



' 7 T/ \ , 

duo Jioj) coth — . 

2kBT 



(95) 



The total transmission coefficient, eq. 



T(e') = ^ rdt [e-^*-^I-^(^'-)]* + 
nV Jo 



is therefore in the broad resonance limit given by 



c.c. 



(96) 



hT Jo 

We note that the polaron shift has canceled out, reflecting that the escape rate out of the 
central site is so fast that any real part environmental self-energy dressing effect is absent. 
The integral in eq. (p6D can now be performed and we obtain 



T(e' 



2r;r 

1 - $ 



r[r-^(e'-e. 



+ c.c. >, 



(97) 



where ^{z) is the probability integral, the error function. We note, that in the above case 
the resonant line shape is no longer of the Lorentz type. The energy dependence in the 
broad resonance limit of the transmission coefficient in the ohmic case, s = 1, is illustrated 
in figure ^. 

b. Strong coupling limit. In the strong coupling limit, rj > 1, we only need to consider 
the short time limit, uJct < 1, because at larger times the influence function is exponentially 
small due to the short range of the correlation part, Z^{t). Therefore, the correlation part 
of the influence function is for the present case the same as in eq.(p^), and similarly the 
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retarded part of the influence function is specified by eq.(0). The expression for the total 
transmission coefficient is therefore the same as the one given in eq.(0). However, the 
vafidity condition here is the strong couphng criteria, 1] > 1, and there is no requirement on 
the escape rate F. We note again, that in the short time hmit approximation the retarded 
part of the influence function has no effect, and the correlation part of the influence function 
is for the present case the relevant one. In flgure ^ the strong coupling approximation for 
the transmission coefficient is plotted versus energy for the case of an ohmic bath. We get 
in the strong coupling approximation for typical parameter values, t] = 10, huc = T and 
2kBT = r, a 5% too low transmission maximum which is displaced to a slightly higher 
energy as compared to the exact result. 

c. High temperature limit. When the temperature is high enough, ksT > (1 + l/r])hiJc, 
we can again use a short time approximation 

Z^(t)=exp|-^t^} (98) 

where now 

,'(r) = H%In.£M. (99) 

h Jo oj 

which again controls the retarded part of the influence function to be 

Z^(t) = expj-z^tj . (100) 
Similar to the two previous cases, the expression for the total transmission coefficient is 



+ C.C. , (101) 



which has the same form as eq. (p7D , except that k, is replaced by k' . The high temperature 
approximation is compared to the exact result in flgure |^, for the parameter values rj = 1, 
huc = r and 2kBT = lOF. We note that the high temperature approximation gives a slightly 
lower transmission maximum than the exact calculation yields. 
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d. Weak coupling and low temperature limit. When the couphng is weak, rj < 1, and 
the temperature is low, fc^T < huc/rj, the situation is different from the previous cases. The 
weak couphng condition, r] < 1, forces the part of the influence function, Z^, to equal unity 
at all times, Z^{t) ~ 1. We also note that there is no contribution from the correlation part 
in the short time limit due to the weak coupling condition. We need therefore to consider 
the long time behavior of the correlation part. The reason for the condition ksT < hcUc/v 
is to avoid the high temperature regime where only short times give a contribution. 

The argument of the exponential in the expression for the correlation part of the in- 
flunence function 

poo J{u}\ tlUJ 

Is{t)= dw^j^coth— —{1- COS ujt) (102) 
Jo oj IkbJ- 

approaches in the long time limit, t > max(^, jj^)^ the expression 

I^{t) = v^^t{u}cty-' / dxx'-^{l - cosx), (103) 
n Jo 

for the exponent region < s < 2. If s > 2, then /^(t) approaches a constant in the long 
time limit, and there is therefore no contribution from the correlation part in this case. To 
be more specific, we perform the integration for the ohmic case, s = 1. In the long time 
limit, Ii{t) = rj^^^^t, and the ohmic correlation part becomes 

ZA-«,=exp{-^(). (104) 

Therefore, the total transmission coefficient is 



or.r roc 



. 1 nkBT 
exp \ - -r,-^t 



~[r - i{e' - Q]t^ + c.c. 



(105) 



The polaronic shift is no longer canceled as ~ 1. Performing the integration we find for 
the total transmission coefficient 
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where 

r., = r + ^. (107) 

We note, that the Lorentzian character of the resonant transmission coefficient is preserved 
in this hmit, but the width and the height of the resonance peak are different compared to 
the case where the bath is absent. Clearly, at high enough temperatures, T > the 
resonance peak will be strongly reduced. Therefore, although the coupling is weak, the result 
is highly non-perturbative. The weak coupling approximation in the temperature range 
^ < ksT < ^ is compared to the exact result in figure |^, for parameter values r] = 0.1, 
huc = lOOr and 2kBT = lOOF. For this parameter choice the resonance peak is broadened 
approximatively five times compared to the case where the environment is absent. Even 
though in the case considered the coupling constant is not that small and the temperature 
is not excessively low, the deviation of the approximately calculated transmission maximum 
compared to the exact value is no more than 10%. 

We notice that the zeroth moment of the transmission coefficient, 

/ de'T{e') = —^, (108) 

is preserved exactly in all the approximations considered above. We also notice that the 
ffist moment of the transmission coefficient, 

^jy,Ti.')=.. (109) 

is preserved exactly by all the approximations except for the weak coupling and low tem- 
perature approximation where the transmission coefficient is symmetric around ec, which is 
anyhow close to ec because the coupling is weak. 



B. The fluctuating level model 

When discussing specific properties of a physical system, it is often possible to neglect 
the quantum nature of the environment and represent the effect of the environment by a col- 
lective classical variable, which is appropriate when quantum fluctuations can be neglected. 
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A well known example is the excellent account of the atomic spectra obtained by disregard- 
ing the quantum fluctuations of the electromagnetic environment, except for cases where 
degeneracies are only lifted by radiative corrections. A counter example where the effect of 
the environment is of pure quantum nature is of course as easily recalled, that of stimulated 
emission. 

In the case where the degree of freedom X represents environmental degrees of freedom 
collectively in the form of an external classical potential, we can still use the presented 
method for calculating the influence. In this case the quantity X just corresponds to a 
classical potential, and is not an operator but just a c number. We shall in this section 
discuss the case where the potential X is a fluctuating quantity, the fluctuating level model. 

We assume further that the fluctuations are Gaussian, and therefore characterized by 
the lowest order correlations 

(X(t)) = c 

{[X{t) - (X(t))][X(t') - (X(t'))]) = K{t - t'). (110) 

The influence of the now classical environmental degree of freedom X is given by the gener- 
ating functional expression, eq. (^4[) , however now the brackets simply denotes the Gaussian 
average with respect to the fluctuating quantity X. We therefore have for the generating 
functional, which in this case is just the probability theory contour characteristic functional, 

^[/] = exp{-^^rfr/(r)(X(r)) 

/ rfr ^ rfr' f{r)K{r - r')f{r') } . (Ill) 

The reason we in the previous section made the distinction in the effects of the environment 
as expressed in the spht 

Z = Z^Z^, (112) 

where we distinguish between the retarded and correlation contributions, is that they rep- 
resent two distinct influences of the environment. The term Z^ represents the systematic 
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friction type and energy renormalization influence and the fluctuating part, including 
thermal as well as quantum fluctuations. For the thermal environment discussed in the 
previous subsection, the two types of influence were not independent, but related through 
the fluctuation-dissipation theorem. 

We now compare the thermal quantum environment with the classical stochastic envi- 
ronment introduced above. It follows directly in the presented real-time approach that the 
fluctuating level model is in one-to-one correspondence with the fluctuational aspect of the 
thermal case with the prescription for the correlation part 

D^{t, t') ^ -2iK{t - t'), (113) 

and substitution of one for Z^, Z^ H-i> 1. In addition, we see that the polaronic shift 
corresponds to the average displacement of the environment 

" " -c, (114) 

obtained from eq. (ffGl) . 

The transmission coefficient for the stochastic environment is obtained from 

zLF F, roo POO 

T{e') = ^^ / dr' / rfrZ^(r,r') 
n Jo Jo 

exp |^(e' - + iT)t - '-{e' - e, - zr)r j (115) 



where 



( 9 roo roo > 

Z^(r, r') = exp --^ / dt dt' f_{t)K{t - t')f_{t') , (116) 
L ri J ~oo J ~oo J 



and 



-ie{T-t)-9i-T' + T-t))] (117) 

and the correlation part Z^ now depends on the explicit form of the correlator K. For 
example, the high temperature form of the correlation function for a thermal ohmic 
environment, J{uj) = rjuj, corresponds to a delta-correlator 
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K(t-t')^'J^S(t-t'), (118) 

the white noise case. The fluctuating level model thus only captures the correlation part 
of the thermal environment and neglects the systematic retarded part, except for its average 
influence. 

For the fluctuating level model we therefore get for the influence function in the Ohmic 

case 

Z^(r, r') = exp ^^-l.n,^kBT\T - r'|| , (119) 



and for the transmission coefficient 



Teff 4r,r. 



Comparing the result to the uncoupled case, we notice, besides the average energy shift and 
a reduction of the peak,Ei that the resonance width T of eg. (PI) is broadened according to 

r ^ Teff = r + ^vksT. (121) 

We note that the transmission coefficient in this case is the same as the one obtained for 
the thermal case in the weak coupling and low temperature limit. 



C. The external field case 

Just as in the previous section we shall here consider the case where the environment 
can be described classically. In addition we shall assume, that in contrast to the previous 
section we know not only the probability distribution of the potential, but in fact the actual 
potential. In the following we shall therefore investigate the situation where we are able 
to couple the electron to an external field at the central site. In the present case of linear 
coupling it is sufficient to consider the case where the central site energy level changes 
harmonically in time. The external potential is therefore given by 

X{t) = Xocosujot. (122) 
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Such a situation could for instance be realized in the case where the Hamiltonian represents 
a small metallic grain with a single active level whose energy can be changed by an external 
electric potential, and the grain being coupled to large metallic electrodes. 

In the case where the influence of the environment is represented by the external potential 
given by eq.( |122D we obtain for the influence function 

= e"^^" /^a!t(/i(t)-/2(i))cosa;ot 



= exp < - — (smcjoti — smci;ot:2 

+ sincjo^s — sinco'o^4)|- (123) 

This corresponds to the fluctuating level model for the case where the potential X{t) is 
known with certainty to be given by the expression in eg . (11221) . 

For the transition probability per unit time we then obtain the expression 

Wp'^i^p.r = ^''^ .4 ^ ' / dr" / dr' / dr 

h J-oo Jo Jo 

exp |-^(ep,r - ep',;)r" 
-^(cp',' - - «r)r' + ^{(^p',i - + »ir)r| 
lim ^ [ dtiexpl^ — (3fJe/^° , ^„ sin uoti 



1 /■* .. (^Xo 
+53m/,^°,_,„coso;oti)}, (124) 
where we have introduced the function 

j!UJo^ — ^ _ gil^OT _j_ ^iLO0(T+T") _ ^iwo{T+T" -t') (125) 



We now calculate the integral over ti in eq.( |124D by expanding the exponential before 



performing the integration, and therefore assume that ujq is non-zero. The case of zero Uq 
corresponds to a constant external potential, X{t) = Xq, which only yields a trivial shift of 
the resonant energy. The transition probability per unit time can now be expressed as an 
integral involving the zeroth order Bessel function, Jq.^ 
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2P^P^ /-oo roo POO 

T(e, e') = 3^ / dr" / dr' / dr 



Tch"^ J~oo Jo Jo 
^°(^'^-°'-"') 

exp{-l(e-e')r"-l(e'-e,-^r)r' 

+ ^(e'-ec + ^r)r|. (126) 

If we had assumed that the external potential was a sine function, X{t) = XosinuJot, we 
would have obtained the same transmission coefficient for non-zero ujq, but we would of 
course not have obtained a shift of the resonant energy in the case of ujq = 0. The total 
transmission coefficient is 

exp |-^(e' - ee - iT)t' + ^(e' - + zr)r| . 

(127) 

Noting that |/^°, ^,,^qP = 4sin^ ^ uio(t-t ) ^ ^ ^^^^ using the summation formula for the Bessel 
function JqB we obtain the total transmission coefficient0 

T(e')=4r,r, Y: n t\2^r2 - (128) 

The total transmission coefficient has Lorentzian side bands at all harmonics of Uujq with 
a relative weight determined by the Bessel functions, and maximal peaks in the spectrum 
at n ~ iXo/ZicJo- The transmission coefficient is plotted for the parameters F = 0.2^ciJo 
and Xq = Ahujo in figure ^. For this parameter choice the first sideband is almost missing 
because is close to the first zero of the first order Bessel function. 

D. Coherent state environment 

In the following we shall investigate the model for the environment in the coherent state. 
Such a model may describe transport through a quantum dot coupled to an environment of 
coherent phonons.0 
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The coherent state has the following representation in terms of the vacuum state, |0), of 
the environment 

|0)=iV;^exp{#S}|O), (129) 

where 

iV^ = exp{|0|2} (130) 

is the normalization factor. 

The state of the environment is then a pure state and the statistical operator for the 
environment reduces to the projection operator 

p, = P^ = (131) 

To calculate the influence of the environment in the case of a coherent state we therefore 
need to evaluate the following expression for the influence function 

Z{UMMM) = Z[f] 

= iV^-i(O|exp{0*6o} 

(Tce-^/c'^-^"(-)^W)exp{06t}|o). 

(132) 

The bracket in the definition of the generating functional corresponds therefore in this 
case to taking the expectation value with respect to the coherent state 

Z[f] = {Tce-rJ.^^n-)^-)) = {(f)\Tce-Uj-f(r)nr)\^^ 

= Ar-i(0|e'^*''« Tce-^/c'^-/W^W e-^^^lO). (133) 



The matrix element appearing in eq.( |133[ ) is most easily calculated by introducing the 
single mode generating functional 

Z[f,f*] = (Tce"^-^='^^^°^-^*^^^^°^^^"^-^^^^^o^^^^). (134) 
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We then notice, that we can rewrite the generating functional of interest as 

Z[f] = {<j,\Tce-^nL^rnr)xir)^^^ 

= Z[fJ*] 



(135) 



provided we substitute into eq.(| 



/(r)= /(r)+z;if 5"(r) 



(136) 



and 



nT)= f{T)-zh^6\r) 

^0 



(137) 



where 5"*^'^ (r) is a delta function on the upper (lower) part of the contour, and vanishes on 
the lower (upper) part. 

The single mode generating functional involves a Gaussian average and is given by the 
quadratic form 



Z[f,r] = e-^Jc-^Jc 



Sc^^ /prfr'A2r(r)iJ(r,T')/(r') 



(138) 



where B is the single mode Green's function 



B{r,T') = -i{Q\Tc%{T)hl{T'm 



(139) 



the zero temperature limit of the previously introduced bath propagator D{t,t') for the 
considered mode. If we therefore insert the proper "force" according to the prescription 
eqs.( p.36| , [137| ) we obtain the influence function for an environment in a coherent state 



ZT=o(t4, ts, t2, ti 

Ao 



exp 



2i 



(140) 
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Here Zt=o denotes the influence function for the thermal case with the temperature set 
equal to zero, that is, the ground state case. 

Exploiting the observations already made for the external fleld case we obtain the fol- 
lowing expression for the transmission coefficient 



T(e,6') = ^/ dr" I dr' I dr ZT=o{r,r' ,t"' 
ttTi J-oo Jo Jo 

T pAo|0||..o A 

exp {-^(e - e')^" - ^(e' - - «r)r' 

+^(6'-e, + zr)r}, (141) 

where the function /^°, ^„ is the same as the function introduced in eq.( |125| ). 

The coherent state case shares features with both the thermal and the external fleld 
cases. The coherent state influence function consists of a factor equal to the zero temperature 
Einstein model influence function, which describes the systematic influence and quantum 
zero point fluctuations, and in addition a factor identical to the external fleld influence 
function. 

For the total transmission coefficient we get 



(l - e-^'^o^"-"')) 



Xo:' 

hoJn 



exp |-^(e' - - iT)r' + '-{t' - + iT)t^ . (142) 

Performing a calculation similar to the one for the oscillatory level model we get a simple 
formula for the total transmission coefficient 

r(e') = 4r,r. £ . , . (i43) 



where 
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In the limit = 0, we of course recover the zero temperature Einstein model result, eg. (pO|) . 
In the limit of Aq ^ 0, ^ oo, and Ao|0| = constant, we recover the classical oscillatory 
level model result, eq.( |128| ). This cross-over behavior is illustrated in figure |^. 

VII. SUPPRESSION OF QUANTUM INTERFERENCE 

In the following section we investigate how the different environments influence the phase 
coherence in a quantum interference setup. We envisage the situation where two resonant 
levels coupled to environments are placed in parallel and transport can take place through 
eitherlll. A physical realization could be double barrier structures situated on the two arms 
of a ring coupled to two reservoirs. The Hamiltonian still has the form 

H = H, + H, + H,, (145) 

but now the sample Hamiltonian corresponds to two levels 

Hs= ^cala,. (146) 

c=l,2 

The tunneling can take place through either level with different couplings 

Ht=Yl {W«I«P,<x + h.c.}. (147) 

p,(T,C 

The interaction takes the form 

H,= ala.X,, (148) 

c=l,2 

where the environment operator depends on the level through the coupling constant 

Xc = T.Kc{bl, + bc.,,}, (149) 

a 

and we assume a situation where the double barrier structures are coupled to separate 
environments 

H,= ^^JblA,c + l]. (150) 

c=1,2;q 
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In order to have an external parameter to vary we envisage an Aharonov-Bohm type 
situation by piercing the ring with a magnetic flux, $, so that the propagators change 
according to propagation around the different arms of the ring according to 

^ e5*/*oGf , ^ e'^^/^^G^. (151) 

The transmission probability given by eg . ([Til) now consists of transmission through either 
arm and to accommodate this two-level situation eq. ([l^) is changed into 



Pp',l^pr{t) - -4 ^P',',ciyp*r,c2^P,'-,C3^p',i,C4 

Ci,C2,C3,C4 



I dti I dt2 / dt3 / dt4, 
JO Jo Jo 

exp |-ep,^(t2 - ts) + ^ep',;(^4 - h] 
{GirAU,h)GUh,h)). (152) 



where 



ih, ti) = -tOih - ti){o\[d,, it2), dl (ti)] |0) (153) 



and 



^t, 03(^4, h) — G^ {ts^ti 



= i9{h -U){0\[d,,iU),dlih)]\0). (154) 

In the following we shall neglect all terms except those for which C4 = C3 and C2 = Ci. This 
is justified if the two resonant levels have an energy difference larger than the width of the 
levels. Here we simply implement it corresponding to propagation taking place through 
either arm. 

The transmission probability given by eq.(|ll]) then consists of the transmission proba- 
bilities for transmission through either arm and an interference contribution: 

Pp',l^pAt) = E Pcip, P', t) + P,nt{p, P', t), (155) 

c=l,2 
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where 



Pc{p, P',t) 



i,r,c 



I dti / dt2 / dts / dt4 
Jo Jo Jo 

exp |-ep,r(t2 - h) + ^ep',K^4 - ^i' 

(t4, t3)Gf (t2, h)Z,{U, t2, h) (156) 



and the interference term is given by 



t rt rt rt 

dti / dt2 / (its / dti 
Jo Jo Jo 

exp |^ep,,.(t2 - ta) + ^ep/,K^4 - 

G^^(^4, t3)G'f (t2, ^l) 

Zi(0, 0, t2, tl)^2(t4, ^3, 0, 0) + C.C., (157) 

We now have an influence function for each of the two different environments, Zi and Z2, 
corresponding to the two different arms of the interferometer. For the transmission coefficient 
we now have according to eqs. (|l8| , [T9|) 

T(e, e') = hm ) Pp',i~^pAtW - ^p'M^ - ep,r-) 

P,P' 

= Tc(e,e')+T,„i(e,e'). (158) 

c=l,2 

The calculation of the total transmission coefficient, specified in eq.(0), in various limits 
is similar to what has been calculated in section 0. For example, in the absence of an 
environment we get for the total transmission coefficient 



T(e') = Ti(e') + T2(e') + 2v'Ti(e')T2(eO 

|ri2r| |r2iH 



e'-ei 



T2 \ 

+ arctan h arg ri2rr2i; , (159) 

e — £2 / 
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where = 27r$/$o is the relative phase difference due to the external flux, and Ti and T2 
are the total transmission coefficients for the individual arms, and specified in the previous 
sections. The decay rates, whose energy dependence can be neglected in the wide band 
limit, are defined as 

Te = T,r + Tel (160) 

and 

rcic2a(e) = TT ^ ypac^^P<rc2^{^ - ^pa) ■ (161) 
P 

To illustrate the special features of the suppression of phase coherence in the present model, 
we discuss the cases of thermal bath, fluctuating level, classical field, and coherent state 
environment below. 



A. Thermal bath and fluctuating level cases 

The transmission coefficient through either arm we obtained previously. For the inter- 
ference term, we have the expression 

2 

Tinti^, e') = e**^ ri2rr2ii 
nn' 

/oo roo too 

dr" / dr' / dT 
-oo JO Jo 

Z^{t, 0,0)Z2(0,r',0) 
exp{-i(6-6')r"} 
exp {-^(e' - - «r2)r'| 

exp|+^(e'-ei+iri)T| + c.c., (162) 
where the influence function corresponding to site c enters as 



Zc{t, 0, 0) = exp <^ -- duj 



a;2 



■ sm ujT 



exp < — 7 / duj — coth 



4 Jo uj"^ 2kBT 
x(l -costur) I, (163) 
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and the central site energies for each arm are shifted downwards according to 

~ec = ec - ^, (164) 
with the negative polaronic energy shift 

X^= Tduj^, (165) 

Jo UJ 

being site dependent through the spectral function 

Jciu;) = ^J2Xljiu-u^). (166) 

In the course of the derivation we have noted that 

Ze(r,0,0) = Z:(0,r,0), (167) 

which follows from eq.(|l^). The interference contribution to the transmission coefficient is 
seen to have the form 

Tint{e,e') 
4 



2 e"f' TurT2iiSie - e') 



J^^dr' Z^iO, r', 0) exp |-^(e' - e~2 - ^r2)r'| 
J^^dr Z*{0, r, 0) exp |^(e' - e\ + tT,)T^ 

+C.C. (168) 

We notice, that the interference term vanishes unless the energy of the particle is conserved. 
Any energy exchange with the thermal environment thus destroys the interference. In the 
present model interference is thus hypersensitive to the presence of the thermal environment, 
and in a dramatic fashion displays the equivalence of dissipation and loss of phase coherence. 

We can also calculate the transmission coefficient for the case of fluctuating levels. We 
assume that each site, c = 1 and c = 2, is coupled to separate classical fluctuating environ- 
ments, represented by the variables Xc{t), c = 1,2, respectively. Analogously to the previous 
thermal bath case there are contributions to the transmission coefficient from transmission 
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through either arm, and there is an interference contribution specified by eq.( |168|) , except 
for the thermal infiuence function being replaced by 

Ze(r, 0, 0) = Z,(0, r, 0) 

= exp I — \ rdt' Kcit -t')\ , (169) 
I Ah Jo Jo J 

where 

K,{t,t') = - (X,(t)))(X,(t') - (X,(t')))). (170) 

As in the thermal case, the interference contribution to the transmission coefficient is com- 
pletely suppressed by dissipation. 

B. External field and coherent state cases 

In the case where the energy levels change harmonically in time with the same frequency 

X^{t) = X^cosujot. (171) 
but with different coupling strength, we obtain for the interference term 
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where we have introduced the function 
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_^lll(^giu;o(T+T") _ ^iujo{T+T"-T')y {^7^) 
TvjJq 

Upon performing the integral over ti we obtain 

Tinti^, e ) = -^^I2r^21l 

exp{-l(e-e')r"-^(e'-e2-ir2)r'} 
exp|^(e' - ei + iri)T| 

+C.C. (174) 

Change in energy of the particle due to interaction with the external field does not destroy 
the interference in accordance with the notion that such a non-dissipativc process does not 
lead to suppression of phase coherence when time reversal symmetry is unbroken. 

Carrying out a similar calculation for the coherent state case, we find a similar qualitative 
behavior for the interference term. Exchange of one quantum does not change the coherent 
state and interference is not destroyed. 



VIII. SUMMARY AND CONCLUSION 

We have developed a functional method to study the influence of various environments 
on quantum tunneling, and shown that the effect of the environment on the transmission 
probability, specified by the influence function, is described by the value of a certain gener- 
ating functional. The state of the environment only occurs explicitly in the expression for 
the influence function through the appearance of its statistical operator, thereby allowing a 
unified discussion of the influence of environments on phase coherence and thereby on the 
transmission properties. An advantage of the approach is that different environments are 
treated on an equal footing, thereby simplifying comparison of similarities and differences. 
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We have calculated the influence function for various environments and parameter 
regimes using the developed non-equilibrium generating functional technique. In order to 
obtain analytical results we have concentrated on simple resonant tunneling, and invoked the 
wide band approximation. In the thermal case we obtained analytical results for the broad 
resonance limit, the strong coupling limit, the high temperature limit and the weak cou- 
pling and low temperature limit. The short time approximation used in the broad resonance, 
strong coupling and high temperature limits, and the long time approximation invoked in 
the weak coupling and low temperature limit were numerically demonstrated to be good 
approximations. For the thermal case we discussed the closed expression obtainable for the 
Einstein model. 

We studied the effect of a classical fluctuating environment on the tunneling process, 
the fluctuating level model. The environment is here described by a Gaussian distributed 
stochastic variable. The fluctuating level model was shown to only capture the correlation 
part of the influence of the thermal environment, as the fluctuating level model neglects 
the systematic part except for its average influence. We flnd in this case that the total 
transmission coefficient is a Lorentzian with its width temperature broadened. 

We studied the effect of a harmonically varying external field, and found the heights and 
widths of the Lorentzian sidebands at all the harmonics of the classical force. 

For the coherent state case we also obtain a closed expression for the infiuence function, 
and we showed that the coherent state case is an intermediate situation sharing features of 
both the zero temperature Einstein case and the external field case. 

Finally we studied the suppression of quantum interference by considering a simple model 
of two resonant levels situated on the two arms of a ring connected to two external reser- 
voirs. The quantum interference was monitored by having the ring enclosing a magnetic 
fiux. When calculating the transmission probability we obtained contributions correspond- 
ing to tunneling through either arm and an interference contribution. We found that the 
interference term vanishes in the thermal bath and fiuctuating level cases even for the slight- 
est energy exchange between system and environment. The model thus represents a case 
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where dissipation completely destroys quantum interference. In the case of a non-dissipative 
environment, e.g., a classical force, there is no loss of phase coherence, and there is therefore 
always an interference contribution, even if energy is not conserved in the transition. 

In conclusion, we have demonstrated the effiency of the developed non-equilibrium gener- 
ating functional technique for evaluating the environmental influence on coherent quantum 
processess. 
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FIGURES 



FIG. 1. The Feynman diagrammatic representation of the decoupling of particle and environ- 
ment dynamics for the amplitude correlator {G^{t4, tz)G^{t2-, t\)) in the wide band approximation. 
The double lines denote the full central site Green's functions while the single lines denote the cen- 
tral site Green's functions in the absence of coupling to the environment. The curly lines denote 
the correlator of the environment operator. 



FIG. 2. Closed time path C extending from — oo to -|-oo along Ci and back again along C2. 

FIG. 3. The total transmission coefficient (in units of ATiTj./T'^) for the Einstein environment 
case. The parameters are V = 0.2hujo, Aq = and ksT = TvjJq. 

FIG. 4. Transmission coefficient (in units of 4r;rj./r^) for the ohmic environment case in the 
broad resonance limit, T > hiOc- In the figure the broad resonance approximation (solid curve) is 
compared to the exact result (dots). The parameter choice is: rj = 1, 2kBT = T and fiLOc = O.IF, 
which yields ^^^(r) = O.IOIF^. 

FIG. 5. Transmission coefficient (in units of 4TiTr/r'^) for the ohmic environment case in the 
strong coupling limit, ij > 1. In the figure the strong coupling approximation (solid curve) is 
compared to the exact result (dots). The parameter choice is: rj = 10, HlOc = F, 2kBT = F and 
n^K{T) ^ 14.27 F^. 
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FIG. 6. Transmission coefficient (in units of 4r;rj./r^) for the ohmic environment case in the 
high temperature hmit, ksT > (1 + l/ri)hujc- In the figure the high temperature approximation 
(soHd curve) is compared to the exact result (dots). The parameter choice is: rj = 1, fiujc = T, 
2kBT = lOr and h'^K'{T) = 10 T^. 

FIG. 7. Transmission coefficient (in units of 'iTiTj./T'^) for the ohmic environment case in the 
weak couphng and low temperature limit, r] < 1 and ksT < In the figure the weak coupling 
and low temperature approximation (solid curve) is compared to the exact result (dots). The 
parameter choice is: rj = 0.1, huJc = lOOT and IksT = lOOF. 

FIG. 8. Transmission coefficient (in units of 4r;rr/r^) for the external field case, 
X{t) = Xq cos ujQt. The parameter choice is: P = 0.2fiujQ and Xq = 4?iwo- 

FIG. 9. Transmission coefficient (in units of 4TiTr/r'^) for the coherent state environment. In 
diagram (a) we have chosen F = 0.2huJo, Xq = Hlvq and = 1. In this case the transmission curve 
resembles the one obtained for the Einstein model, see figure |3[ In (b) F = 0.2hu!Q, Xq = Hujq and 
|(/>| = 2. In (c) F = 0.2hujQ, Xq = 0.2hujQ and = 10. The shape of the curve in diagram (c) 
resembles the one obtained in the external field case, see figure ^. 
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